%% =========================[comparison]================================
run quasi0StiffSuspension.slx
close all
step=0.2;
time=0:step:10;
timeNumber=length(time);
accAlone=zeros(timeNumber,1);
for i=1:timeNumber 
    for j=1:length(acc.time)
        if acc.time(j)==time(i)
        accAlone(i)=acc.data(j);
        end
    end
end
plot(time,accAlone,'-*','linewidth',1)
hold on
load MCdata
errorbar(time,PCmean,sqrt(PCvar),'-h',...
    'color','black', ...
    'linewidth',1,...                        
    'MarkerSize',6,...                        
    'MarkerEdgeColor','red',...               
    'MarkerFaceColor','white')                   
title('Mean and standard deviation of acceleration over time')
xlabel('Time (s)')
ylabel('Acceleration (m/s^2)')
